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. SUMMARY 


The widely used transonic swept-wing code, FL022, was found to have an error 
in the transformed flow equation in the computational domain. The revised version of 
the code correctly accounted for the non-straight-leading-edge geometry and its 
effect on the pressure distribution. 


INTRODUCTION 


Transonic flow calculations have been made for many years. Among the wide 
variety of numerical procedures that are now available, a series of computer pro- 
grams known as the FLO codes have been developed by Professor Jameson and his 
coworkers. One of them, the widely used swept wing code, FL022 (ref. 1), has been 
known to generate inaccurate results for wings with non-straight-leadlng-edges . 
Recently, the code was used to calculate the flow past an ONERA cranked wing model 
mounted on a wall (ref. 2), and the predictions were found to be very poor in the 
vicinity of the crank. The reason for the poor predictions was identified as an 
error associated with the scaling between the physical space and the computational 
domain. In fact, the FL022 code is dimensional for wings having non-straight- 
leading-edges. 


ANALYSIS 


The governing equations (ref. 1) for steady, Inviscid, compressible isentroplc 
flow in Cartesian coordinates (x,y,z) are the potential-flow equation 

(a^ - u^)(i) + (a^ - v^)(i) + (a^ - w^)((i - 2uvd> - 2uwd> - 2vw(j) =0 (1) 

XX yy zz xy xz yz 

the Bernoulli equation 


j (u^ + v^ + w2) + 


= constant 


and the isentropic relation 


= constant 


Here, (f> is the velocity potential (defined . such that V(j> - (u,v,w) is the velocity 
vector), a is the speed of' sound, y is 1.4 for air, p is the pressure, and p is 
the density. Let the reduced potential, G, be defined through the relation 

, 4 ^ - r. L *• ■ 

• G = <J> - X cos a - y sin a 



The term G satisfies equation (1), where a is the angle of attack of the wing. 
(Far-field flow speed and density are taken as unity.) 


The flow equation (1) can be written in terms of G and solved in a computa- 
tional domain generated through a sequence of mappings from the physical space. 

The first three mappings consist of the first shearing transformation 

X = X - x^(z) 

y = y - yg(z) > ( 2 ) 

Z = Z 


the square-root transformation 

(Xi + iYi)2 = 2(x + iy) 
Zi = z 

and the second shearing transformation 

X = Xi 

Y = Yi - S(Xi,Zi) 

Z = Zi 


(3) 


where the singular line (the Jacobian vanishes) of the square-root transformation is 
represented as x = Xg(z), y = y^iz) and the entire wing surface is described as 
Yi = S(Xi,Zi). 

It is useful to introduce a scaling mapping between the first shearing transfor- 
mation and the square-root mapping 


(4) 



that will account for the exact scaling relation between the physical space and the 
computational domain. The parameters SCAL and SCALZ are real variables that are 
defined in the subroutine COORD of FL022, and set up the stretched coordinates 
(see ref. 1). The square-root transformation becomes 


(Xi + iYi)^ = 2(x + iy) 
Zi = z 


} 


(5) 


The flow equation, in terms of G and using the mappings (2), (3), (4), and (5), 
becomes an equation of the form 


2DG^ + 2EG^^ + 2FGy, 


+ RiG^ + R2G^ = 0 


( 6 ) 
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which is equivalent to equation (9) of the original analysis (ref. 1) if the follow- 
ing notations are introduced; 



C = -(x^o + y^y) 
n = x^y - y^o 
a = -(oSx + y) 

6 = 0- ySjj 

Q = SCAL/SCALZ 


Y = n - CSj, - QS^ 

X = x'Xi~~ + yIXi~~ 
s -ixx ■'s J-xy 

'V = -ylXi-- + x’Xi~~ 

•'s ^xx s -^xy 

A = SxXi~~ + Xi~- 
^xx ^xy 

Z = SxXi~- + Xi~~ 

^ ^xy ^xx 

U = uo + vy + wC 
and V = ua + v6 + wy 

Then the coefficients of equation (6) can be written as 


A = (o^ + y2 + _ \j2 

B = (a2 + g2 + y2)a2 _ y2 

C = Q^(a^ - w^) 

D = (oa + y6 + Cy)3^ " UV 

E = Q(5a2 - Uw) 

F = Q(Ya^ - Vw) 

Rl = -SCAL(a^ - w^) (x''o + y"y) + (x'x + y'l') - Xi~~(u^ - v^) - 2uvXi-~ + 2uwx 

s o s s XX xy 

+ 2vw'i' 


and R 2 = -SCAL(a2 - (x"a + y"6) + U^S + w^Q^S + 2UwQS + (u^ - v^)A + 2uvZ 

S S XX zz xz 

- 2uw(S^x + t) - 2vw(s^'i' - x) - a2[(a2 + y2 + + 2CQS^^] 

+ (a2 - w2)[(y^X - x^f) - S^(x^x + 
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In the original version of the FL022 code, the real variable SCAL was omitted 
in the expression for and R 2 , which is equivalent to setting this variable 

to unity. This did not have any effect for straight leading-edge geometries, but 
caused errors for non-straight-leading-edges, or those with curvature. 


RESULTS 


The flow over the ONERA cranked wing (fig. 1) was calculated at a Mach number 
of 0.85. The reflection-plate model had an NACA 0012 airfoil and was set at an 
angle of attack of zero degrees. Comparisons with the ONERA experimental data and the 
original FL022 results are shown in figures 2-6. The correction improves the corre- 
lation significantly near the leading-edge crank (station #3). 


CONCLUSIONS 


The widely used transonic swept-wing computer code, FL022, was found to have an 
error in the transf ormed-f low equation in the computational domain. The revised 
version of the code correctly accounted for the discontinuous leading-edge geometry 
and its effect on the pressure distribution. 
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ONERA WING ON WALL 


FMACH = 0.85 ALPHA = 0 





Figure 1.- The ONERA cranked wing on wall. 
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STATION = 1 

MACH = 0.8500 AL = 0 

CL = 0 CD = 0.0188 CM = 0 


Figure 2.- Pressure distribution at span station 1. 
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STATION = 2 

MACH = 0.8500 AL = 0 

CL = 0 CD = 0.0118 CM = 0 


Figure 3.- Pressure distribution at span station 2. 


7 


- 1.2 - 
-.8 - 


-.4 - C 



.4 - 


.8 - 


1.2 L 


Figure 4 
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STATION = 3 

MACH = 0.8500 AL = 0 

CL = 0 CD = 0.0047 CM = 0 


- Pressure distribution at span station 3. 
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STATION = 4 

MACH = 0.8500 AL = 0 

CL = 0 CD = -0.0045 CM = 0 


Figure 5.- Pressure distribution at span station 4. 


9 


-1.2 r 


□ 


O 


ORIGINAL 
REVISED 
ONERA DATA 



1.2 



STATION = 5 

MACH = 0.8500 AL = 0 

CL = 0 CD = -0.0222 CM = 0 


Figure 6.- Pressure distribution at span station 5. 
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